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Abstract. We discuss the determination of elemental abundances from high resolution X-ray data. We emphasize 
the need for an accurate determination of the underlying temperature structure and advocate the use of a line 
ratio method which allows us to utilize, first, the strongest lines observed in the X-ray spectra, and second, 
lines that span a rather wide temperature range. We point out the need to use continuous emission measure 
distributions and show via example that modeling in terms of individual temperature components yields errors 
of more than 50%. We stress the need to derive differential emission measure distributions based on physical 
assumptions and considerations. We apply our methods to the Chandra LETGS spectrum of Algol and show that 
nitrogen is considerably enhanced compared to cosmic abundances by a factor of 2 while carbon is depleted by 
at least a factor of 25. Iron, silicon, and magnesium, are all depleted compared to cosmic abundances, while the 
noble gas neon has the relatively highest abundance. 
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1. Introduction 

The determination of elemental abundances in a variety 
of astrophysical objects belongs to the most important 
tasks of observational astronomy and the understand- 
ing of the evolution of chemical elements with cosmic 
time is among the central themes of modern astrophysics. 
Elemental abundances can be measured cither from ab- 
sorption spectra of stellar atmospheres or from an analy- 
sis of the line emission spectrum of nebular emission. In 
both cases the temperature structure of the emitting ob- 
ject must be known before elemental abundances can be 
determined because for a given set of abundances plasma 
temperature is (often) the most important factor control- 
ling the ionization equilibrium and hence the amount of a 
given type of material, say Ovu or Fexvn, in an astro- 
physical object. 

Reliable determinations of chemical abundances are 
carried out from high resolution spectra. While elemental 
abundance determinations of stellar photospheres can also 
be made from a set of suitably chosen filters, abundances 
determined from high-resolution spectra are thought to be 
much more accurate and far less model-dependent. This 
also and specifically applies to X-ray data. The energy 
losses of hot thermal plasmas with temperatures above « 
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10 6 K peak in the X-ray range and therefore the chemi- 
cal composition of such plasmas, which are encountered 
in stellar coronae, in supernova remnants, and clusters of 
galaxies, can in fact only be determined from X-ray data. 

Continuum energy losses dominate the cooling of ther- 
mal plasmas above rs 10 7 K, while the thermal energy 
losses of plasma with temperatures below ~ 10 7 K are 
dominated by a multitude of emission lines. At those tem- 
peratures the strongest coolants are typically (albeit not 
exclusively) the hydrogen- and helium-like ions of the most 
abundant species, i.e., carbon, nitrogen, oxygen, neon, 
magnesium, silicon, sulfur, calcium, argon, and iron. The 
hydrogen- and helium-like ions of the elements heavier 
than sulfur are located below 6 A, and are therefore dif- 
ficult to observe with high spectral resolution. Also, at 
temperatures below « 10 7 K, the dominant stage of ion- 
ization for the heavier atomic species is not yet advanced 
to helium- or hydrogen-like ions. For example, for iron the 
most abundant stage of ionization at 10 MK is boron-like 
Fe XXIII, at a temperature of 3 MK neon-like Fexvn (cf. 
Arnaud and Rothcnflug 1985), and consequently the en- 
ergy losses from iron are dominated by line emission from 
those ions. 

The wealth of emission lines from trace elements in 
the X-ray range demonstrates the potential of abundance 
determinations from such data. Many calculations of the 
total energy output of a hot collisionally ionized plasma 
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have been made (cf. Raymond and Smith 1*11)771 Mewe et 
al. 119851 Dere et al. I1997fl . and the results of such cal- 
culations have been used to interpret broad-band X-ray 
data as available from the Einstein Observatory IPC or 
the ROSAT PSPC and HRI. Such plasma codes have been 
used to infer an energy flux from a given count rate mea- 
surement as well as to model the typically low-resolution 
pulse height spectra of proportional, gas scintillation or 
CCD detectors. For the spectral modeling individual spec- 
tral components are used. While the lowest resolution data 
can be modeled with one or two temperature components 
with solar abundances, higher resolution spectra require 
three or more temperature components with typically non- 
solar abundances. X-ray detectors tend to work very ef- 
ficiently at energies of « 1 keV, where both the effective 
area of the instrument and the plasma emissivity from 
iron is at maximum. Deviations from solar abundances in 
the X-ray spectra of stars were first reported on the ba- 
sis of ASCA CCD spectra of stars like Algol (cf. Antunes 
et al-HHSUl and AR Lac (cf. White et al. lT3ikj) : also the 
even lower resolution PSPC spectra of some active stars 
were found to be better fitted with metal-depleted rather 
than solar-abundance plasma models (for example, Algol, 
cf. Ottmann and Schmitt ITMol and CF Tuc, cf. Schmitt 
et al. I1996[) . During large flares abundance changes were 
inferred on the basis of spectral modeling of the time- 
dependent X-ray spectra. Both in AB Dor (cf. Giidel et 
al. l2UUl[> as well as Algol (cf. Favata and Schmitt the 
iron abundance was found to increase during the early rise 
phase of a flare, and then to decrease back no "normal" 
sub-solar abundance values. 

Abundance determinations of stellar coronae based on 
an analysis of individual emission lines were first carried 
out with data from the spectrometers on board the EUVE 
satellite. The emission line studies based on EUVE data 
followed relatively closely the example of abundance de- 
terminations of the solar corona. Stern et al. I|i595) and 
Schmitt et al. I|199 6 ) found an anomalously low iron abun- 
dance in the EUVE spectra of Algol and CF Tuc, based 
on an analysis of the Fexx, Fexxi, and FeXXll emission 
lines in the XUV range and the observed continuum val- 
ues. Schmitt et al. Ijl996|l coined the term metal abun- 
dance deficiency syndrome (MADS) for this phenomenon, 
which is in contrast to the abundance pattern observed 
in the solar corona, where elements with low first ioniza- 
tion potential (FIP) are found to be enhanced. Drake et 
al. (|1996f> studied the presence and absence of this so- 
called FlP-effect in a small number of nearby stars. Using 
Chandra HETGS data Drake et al. pUtjTfl study the ele- 
mental abundance of the active binary HR 1099 by means 
of a differential emission measure distribution computed 
from a Markov chain, while Brinkman et al. H2()()l(l study 
the same source with the XMM-Newton reflection grating 
spectrometer (RGS) assigning all of the emission measure 
to the temperature corresponding to the peak of the line 
contribution function. Other abundance studies based on 
Chandra or XMM-Newton include Audard et al. fMvty . 
who use a fit approach based on Chebyshev polynomials 



and Huenemoerder et al. Q2001[l . who use a smoothed pos- 
itive emission measure distribution function. Giidel et al. 
(2001) use a fitting approach based on individual temper- 
ature components to derive the elemental abundances in 
YY Gem using again data from the XMM-Newton RGS. 

The purpose of this paper is to apply the ideas devel- 
oped in solar and stellar ultraviolet emission line studies 
to the now available broad band and high spectral reso- 
lution X-ray data. The specific advantage of those data 
coming from the recent generation of X-ray spectrometers 
on board Chandra and XMM-Newton is that they cover 
the resonance lines of the hydrogen- and helium-like ions of 
the elements carbon, nitrogen, oxygen, neon, magnesium, 
and silicon. Resonance lines of the hydrogen- and helium- 
like ions are among the strongest observable lines and can 
be detected also in the weaker sources. Also, the atomic 
physics required for an interpretation of those lines is - 
probably - simpler than that required for the lines in more 
complex ions. We further follow quite strictly the approach 
that the determination of the emission measure distribu- 
tion must occur in an abundance independent fashion, i.e., 
one either uses only lines from a given element (in practice 
only iron can be used) or one uses ratios of lines from the 
same element. The latter approach has the enormous ad- 
vantage that the strongest (rather than the weaker) lines 
in an observed X-ray spectrum can be used for abundance 
determinations, once the overall continuum emission level 
(or possibly that of a well defined atomic species) has been 
fixed. 

The plan of our paper is therefore as follows: We first 
collect the necessary formulae required for the calculation 
of line and continuum fluxes from an isothermal plasma 
at temperature T with specific emphasis on the abun- 
dance dependence of these quantities. We introduce the 
concept of the differential emission measure (DEM); the 
DEM distribution is modeled by an approximation with 
Chebyshev polynomials on the one hand and with the help 
of a Gaussian distribution of magnetic loops on the other 
hand. The abundances computed in this fashion are juxta- 
posed to those computed from a more conventional anal- 
ysis with individual temperature components. 

2. Abundance determination from optically thin 
plasma emission 

2.1. Differential emission measure distribution 

2.1.1. Line emission 

In this section we review the basic physics of coronal line 
formation in as much as relevant for elemental abundance 
determinations. The basic ideas of analysis have been sum- 
marized by Pottasch (|1965|l in a solar context. Consider 
the simplest case of a two-level atom in coronal equilib- 
rium. Coronal equilibrium implies an equilibrium between 
collisional excitation from the lower level I followed by ra- 
diative de-excitation from the excited level u. The emitted 
photon leaves the system so that in essence energy from 
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the thermal pool has been converted into radiation. The 
equilibrium condition then reads as follows: 



i.A 1 . 



(1) 



where n e i, n ion i, and n ionu denote the number densities of 
electrons and ions in lower and upper state, respectively, 
Ciu is the collisional rate coefficient and A u i the Einstein 
for spontaneous radiative deexcitation. The total power 
emitted in the transition ul is given by 



Pal nj /0nu A u iV c 



corona 7 



(2) 



where V corona denotes the total (isothermal) coronal vol- 
ume. Let us further assume that the considered transition 
is produced by a Z-times ionized atom from atomic species 
X; in a low density environment almost all of these ions 
will be in the ground state I, i.e., ni on ,l — n x + - We can 
therefore express ni on ,i in the following way: 



z+ 

_ n£ n x n H 
nx riH n e i 



(3) 



where nx denotes the total number density of all ions of 
atomic species X and tih the hydrogen/proton density. 
The abundance of species X relative to hydrogen is de- 
noted by Ax and is given by 



A x = 



n x 
n H 



(4) 



and in a more or less fully ionized cosmic abundance 
plasma we have 



n p = 0.85 n e i 



(5) 



Note that for all plasmas with "reasonable" composition 
most of the electrons come from hydrogen and helium, so 
that Eq.JHlis virtually independent of the assumed elemen- 
tal abundances. With these definitions the emitted power 
Pal can be written as 



Pal 

— 0.85 n e i Vcorona 



n z+ 

Ciu A} 

n x 



(6) 



The product C; u -^- depends only on temperature; we 
can therefore write 



Pui = A ul (T)EMA x , 



(J) 



where we have defined the volume emission measure EM 
in the usual way through 



EM = ni,V c 



e l * corona 



(8) 



and denote by A u i(T) the so-called line cooling func- 
tion of the transition ul. Eq. illustrates the fundamen- 
tal dilemma of all abundance determinations using op- 
tical thin emission in coronal equilibrium: Cooling func- 
tion, emission measure, and abundance all enter in the 
same multiplicative fashion, and therefore one has to know 
both emission measure and temperature in order to deter- 
mine the abundance. The temperature dependence of the 
collisional excitation coefficients is typically of the form 



Ciu ~ T~ 1 / 2 e~ x ' kT , with \ denoting the line's excitation 
potential, and similar functional dependencies apply to 
the temperature dependence of the fractional ionization 
of a given ion. Consequently, a given line is sensitive over 
a relatively broad range of temperature with a width of 
typically 0.3 dex. 

So far we have assumed isothermality. Stellar coronae are 
very likely not isothermal, and different coronal volume 
elements contributing to a given observed line flux in the 
transition u to I have in general different temperatures. 
Consider an infinitesimal volume element dV in a temper- 
ature range dT, and the differential power dP u i emitted 
from this volume clement. The total power P u i can then 
be computed in the following way: 



P„ 



dP,. 



A x 
A x 



dVn 2 el A ul (T) = 
dT£(T)A u i(T) , 



(9) 



where temperature T was used to replace the integra- 
tion variable V. The quantity £(T) is defined through 



e(T) = n^(T)g 



(10) 



and is known as the so-called differential emission measure 
(DEM) . Obviously we must know the differential emission 
measure £(T) in order to determine the abundance A x - 
Eq. shows, however, that the observed line flux P u i is 
an integral of the product of the line cooling function and 
the differential emission measure distribution. If a number 
of different lines with cooling functions A u i(T) with dif- 
ferent temperature sensitivity are available, one obtains 
an integral equation with the differential emission mea- 
sure distribution £(T) as kernel. The system is a Fredholm 
equation of the first kind. Such integral equations are no- 
toriously difficult to solve and, in particular, need not have 
unique solutions. This has been known to mathematicians 
for a long time, Craig and Brown l|1976|l were the first 
to point out the resulting limitations for our ability to 
extract physical information from spectral analysis in a 
cogent form. 

2.1.2. Continuum emission 

We now consider continuum radiation from a hot 
plasma. Continuum emission comes from (free-free) 
bremsstrahlung, from free-bound radiation, and two 
photon radiation; for temperatures > 5 MK thermal 
bremsstrahlung is the dominant continuum energy loss 
mechanism. All these continuum emission processes orig- 
inate from interactions of either protons or ions with free 
electrons very similar to the generation of line emission. 
One thus expects that the dependence of the continuum 
emission on electron density and temperature is of the 
same functional form as for line emission. Indeed, Mewe 
et al. (|1986|l consider an isothermal plasma with electron 
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number density n e at some temperature T and write the 
specific continuum emissivity as 



dW 



2 5 7re e 



dVdtdX imc 2 



2ir 9 1 he 



(11) 



here h and k denote Planck's and Boltzmann's constants, 
respectively. The quantity g c (X, T) is the velocity-averaged 
Gaunt factor. It can be written as a sum of three compo- 
nents through 



9c = 9ff +9fb + 92 1 , 



(12) 



where gff, gfj,, and represent the contributions from 
free-free, free-bound, and two photon radiation. 

The different constituents of the continuum emissiv- 
ity do, however, have different dependences on elemental 
abundances. The bremsstrahlung component comes pre- 
dominantly from electron-proton collisions with most elec- 
trons due to fully ionized hydrogen and helium atoms. 
Under "reasonable" abundances the number of electrons 
from all heavier elements will be very small. In contrast, 
the other components do depend on trace element abun- 
dances since the free-bound radiation depends on the re- 
combination frequency and the two-photon radiation on 
the number density of meta-stable He-like states, both of 
which are abundance dependent. Fortunately, for a suffi- 
ciently hot plasma electron bremsstrahlung is dominant so 
at least in first order the continuum energy loss is abun- 
dance independent. The specific power dP\ emitted from 
a volume element dV is thus given by 



dPx = 



2 5 7re 6 



2tt 



1 



;<e-w —g c (\,T)dV. 



(13) 



The specific power emitted from a continuous emission 
measure distribution described by a distribution function 
£(T) is then given by 



Pa = 



2 5 7re b 
3mc 2 



^e-^. 9c (A,lW)dT. (14) 



Thus from a mathematical point of view, the contin- 
uum emissivity is again given by an integral of the differ- 
ential emission measure distribution function £(T) with 
yet another kernel. In contrast to line emission, however, 
the continuum emission does, first, not depend sensitively 
on the plasma's metal abundance, and second, its temper- 
ature dependence is very different. At high temperatures, 
where bremsstrahlung dominates, the continuum emission 
is virtually independent of metal abundance since almost 
all electrons in a plasma come from the lightest elements. 
These bremsstrahlung continua are rather flat, the domi- 
nant feature being the thermal cutoff at short wavelengths, 
where the gratings on board Chandra and XMM-Newton 
have small effective areas. Therefore the continuum's ther- 
mal cutoff constrains the highest temperatures existing in 
a corona, but to far a lesser extent the actual run of emis- 
sion measure vs. temperature. 



2.2. Constraints on and modeling of the emission 
measure distribution 

The temperature structure of a stellar atmosphere can be 
derived from the principles of radiative and hydrostatic 
equilibrium, the temperature structure of a magnetically 
confined plasma can be computed from the energy equa- 
tion, if one assumes - for example - a static equilibrium. 
The difficulty in the latter case is, that, first, this tempera- 
ture structure of an individual coronal feature is virtually 
independent on the form of the assumed heating, which 
in essence is unknown, and second, that in the stellar case 
one is very likely looking at the integrated emission of a 
large number of individual features. In other words, this 
integrated emission has to be described by some distri- 
bution function of the physical parameters characterizing 
individual coronal features, and again that distribution 
function is unknown. In consequence we conclude that 
the mathematical form of the differential emission mea- 
sure distribution is a priori unknown. 

Which constraints can nevertheless be imposed on the 
function describing the temperature structure, £(T) ? In 
the following we formulate three conditions and discuss 
their physical and mathematical implications. Clearly, 
from the definition of £ (T) it follows 

€(T)>0. (15) 

This condition looks trivial, but actually represents a 
rather strong constraint. Next, the DEM distributions of 
individual magnetic loops have (integrable) singularities 
at the loop top; if such DEMs are integrated over dis- 
tribution functions, continuous, smooth emission measure 
distributions are found. We therefore assume that £(T) is 
a smooth function of temperature. And finally, we know 
from the many low resolution X-ray data from propor- 
tional counters that the X-ray temperatures cannot ex- 
tend to arbitrarily high values. We therefore assume the 
existence of some maximum temperature T maxi above of 
which no emission measure is present, i.e., 

£(T max ) = . (16) 

In the absence of any plausible physical model we as- 
sume that the differential emission measure £(T) can be 
approximated by a sum of Chebyshev polynomials; wc will 
discuss a physical ansatz in section 5.3. Of course, any 
other system of orthogonal polynomials or other func- 
tions could also be used. Chebyshev polynomials have 
unique normalization properties, and expansions of the 
target function into Chebyshev polynomials are useful as 
long as the expansion coefficients diminish rapidly with 
increasing order. Therefore this approach has also been 
taken, for example by Lemen et al. 1)1989(1 and Stern et 
al. 1)1995(1 . Also, our available X-ray data are only sensi- 
tive to plasma hotter than some temperature T m ^„, since 
plasma with T < T m i n emits outside the band pass. Let 
us then introduce the dimensionless temperature variable 
x defined in the closed interval [0,1] through 
T — T 

- 1 ± mm 



Tn 



T 



(17) 
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Let T n (x) denote the Chebyshev polynomial of order n, 
satisfying the following boundary conditions (n = 0, 1, ...): 

T 2n (0) = (-l) n T 2 „ +1 (0) = T B (1) = 1. (18) 

We seek the differential emission measure distribution 
and approximate it through a sum of M Chebyshev poly- 
nomials via 



A I 



(19) 



i=0 



with coefficients at, i — ... M to be determined from the 
data and the boundary condition 



e(i) = o, 



(20) 



which is equivalent with Eq. 1161 Since we are primarily 
interested in the shape of the emission measure distribu- 
tion (and not in its normalization) we further impose the 
constraint: 



dx £(x) = 1 . 



(21) 



Because of the properties of Chebyshev polynomials 
(cf. Eq. I18fl the boundary and normalization conditions 
Eas. 1201 and 1211 translate into two conditions for the coef- 



ficients Oj, i = . . . M: 



J2 a i = 

f=0 

and 



M .i M 

J^aj / dxTi(x) = J^aj/j = 1 

i=0 J ° i=0 



(22) 



(23) 



with known coefficients 2j. Therefore the number of in- 
dependent coefficients Oj is M — 2 and we can write the 
coefficients clq and a\ as 



M 



and 



i=2 



a = ^(1 - 2 x J^aj - 2 . 

i=2 



(24) 



(25) 



Consider now a specific spectral line denoted by some in- 
dex j with line cooling function Aj(x). The measured flux 
fj of that line is given by 



fj = 47]7 j2 ( T max - T m i n ) A X / dxAj(x)£(x) (26) 



e -*,NB M 1 

-(T max -T rnin ) A x y.a, / dxAj(x)Ti(x) , 

'o 



And 2 



i=0 



where d denotes the distance to the star, Njj the hy- 
drogen column density, and <jj the effective absorption 



cross section per hydrogen atom at the line's energy. With 
the abbreviation 

Lji = [ dxAj(x)Ti(x) (27) 
Jo 

as Chebyshev line contribution coefficients we can write 



fj = 



And 2 



(28) 



i=0 



and one needs to determine M + 1 coefficients <Zj, i — 
. . . M from the data. 

2.3. Differential emission measure reconstruction from 
lines 

A given line flux /,• depends both on temperature and 
abundance (cf. Ea. l28f) . but the ratio of two emission lines 
from an ion of the same atomic species is clearly indepen- 
dent of the specific abundance of the chosen element. In 
order to distinguish between abundance and temperature 
effects on the differential emission measure distribution 
one should therefore work only with line ratios from the 
same elements (or only with lines from the same element, 
which is feasible only for iron). If we let the index z de- 
note an emission line (or possibly a sum of emission lines 
of a given element) in the numerator and analogously the 
index n the line in the denominator, we can write for the 
expected (abundance-independent) line ratio p zn 

fz_ _ J2i=o 

In 



P 



X^ M T 



(29) 



Note that the idea of using line ratios is not new; for ex- 
ample, Mcintosh ( 2000J has used line ratios in an attempt 
to determine the emission measure distribution of features 
in the solar corona, and Fludra and Schmelz l|1995(l have 
used a similar approach in their analysis of X-ray flare 
spectra obtained with the BCS. In our context the great 
advantage of using line ratios lies in the fact that in this 
fashion the helium-like and hydrogen-like ions of many 
atomic species can be used. These emission lines tend to 
be strongest and the most easily detectable ones in any 
stellar corona. For the lighter elements like carbon, ni- 
trogen, and oxygen these are in fact the only detectable 
lines in the X-ray range, and furthermore, atomic physics 
uncertainties should be smallest for such lines. 

We use the double index (z, n) to denote a specific line 
ratio and the expressions Z Z i and N n i to denote the line 
contribution coefficients from the numerator and denomi- 
nator lines entering the ratio (z,n), respectively. Given a 
set N of measured line ratios r Z) „ and errors a z . n , all of 
which are derived from emission lines of the same atomic 
species, we can determine the differential emission mea- 
sure distribution by minimizing the test statistic x 2 de- 
fined as 

2 



E 

(z,n) 



(30) 
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with respect to the expansion coefficients dj for a given 
set of measured line ratios and an arbitrary number M of 
expansion coefficients. Some care is required since these 
coefficients are not independent; we use the boundary con- 
ditions Eas. 1241 and 1251 to express the coefficients ao and 
di in terms of the remaining coefficients aj, i = 2 . . . M. 
We further note that multiplying any set of coefficients 
by some common scale factor T leaves all line ratios P[ z . n ) 
invariant. As a consequence we can express the test statis- 
tic x 2 with the known parameters C z % for the numerator 
and D n i for the denominator in terms of the independent 
coefficients a,- , i = 2 . . . M as 



£ 

(z,n) 



(31) 



For the latter expression derivatives can be computed an- 
alytically and minimization is straightforward. In order to 
ensure that the positivity constraint Ea. ll5l is satisfied, we 
introduce a penalty function with large positive values for 
X 2 once £ (x) becomes negative at any point in the interval 
are [0,1]. 

2.4. Differential emission measure reconstruction from 
continuum 

The reconstruction of a differential emission measure dis- 
tribution from a set of continuum measurements Cj with 
associated errors dj follows in analogy to the line treat- 
ment, except that obviously the appropriate continuum 
contribution functions must be chosen. In general the tem- 
perature sensitivity of the X-ray continuum is much less 
pronounced than that of individual emission lines. The 
most striking feature of the bremsstrahlung spectra is the 
exponential cutoff at short wavelengths resulting from the 
exponential decay in the number of available electrons at 
some given plasma temperature T, while at lower tem- 
peratures recombination and two photon continuum can 
dominate at specific wavelength bands (see discussion by 
Mewe et al. 119861) . We also carried out test calculations of 
continua with "reasonable" differential emission measure 
distributions and found the spectral shape of these con- 
tinua quite insensitive to variations in the parameters of 
the differential emission measure distribution. In our fits 
we therefore used only a small number (3) of continuum 
bands. 

3. Which lines to use ? 

An inspection of Eq. [5] shows that the temperature de- 
pendence of the power emitted in a given spectral line is 
determined by two factors, i.e., first, the ionization equi- 
librium of the atomic species considered and second, by 
the electron excitation rate. The ionization equilibrium 
will in general be such that the fractional ionization of a 
given stage of ionization, say, O VII, peaks at some tem- 
perature, and at lower (higher) temperatures the predom- 



inant ionization will be lower (higher) and the emitted 
line flux will correspondingly change. As to the electron 
excitation rate, in general for excitation of a given line 
a certain threshold energy is required, and all electrons 
above this threshold value will be able to perform atomic 
excitations. Therefore, the excitation rates will in gen- 
eral increase with increasing temperature but eventually 
level out. Therefore, the emissivity of a given line of some 
species, say Ovn, will peak at some temperature T max . 
Note that the line of the hydrogen-like ions are broader 
(with high temperature tails) than the corresponding lines 
from the helium-like ions, an effect caused by the ioniza- 
tion equilibrium. In order to illustrate this behavior we 
plot in Fig. ^ the line emissivity (per unit emission mea- 
sure) as a function of temperature for the Ly a and He-like 
resonance lines (i.e., the transition 1 P\ - 1 Sq) for the el- 
ements carbon, nitrogen, oxygen, neon, magnesium, and 
silicon according to the calculation by Mewe et al. (1985). 
As obvious from Fig. ^ these lines very nicely sample the 
temperature range between w 10 6 and w 2x 10 7 K. These 
lines belong to the strongest lines of these elements in this 
temperature range and are easily observable in Chandra 
LETGS spectra. 

If one now considers line ratios from the same chem- 
ical element but from two adjacent stages of ionization, 
say from O vn and O VIII, the lower ionization stage line 
will dominate at lower temperature and vice versa, and 
the line ratio will monotonically increase with tempera- 
ture. In Fig. [2] we show the ratio p of the energy flux in 
the Ly a line divided by the He-like r line as a function of 
temperature T. As is clear from Fig. the line ratios p 
do indeed increase monotonically with temperature. The 
temperatures at which p is unity for a given atomic species 
increase with increasing atomic mass (and nuclear charge) 
reflecting the fact that more and more energy is required 
to establish, say, the He-like stage of ionization, p increases 
by « 100 for a temperature increase of « 0.5 dex, thus p is 
rather temperature sensitive. Therefore a given measured 
line ratio can be (uniquely) converted into a temperature, 
however, different line ratios will in general result in differ- 
ent temperatures. These temperatures must not be inter- 
preted as "isothermal" temperatures, but rather as "mo- 
ment" temperatures since they depend on the differential 
emission measure distribution (a stellar property) and the 
lines' emissivity functions (an atomic property). 

How do these curves compare to observations? We in- 
cluded the measured data points for the stars Algol (high 
activity), e Eri and Capella (intermediate activity), and 
Procyon (low activity) taken from Ness et al. 2002b who 
gathered together line ratios of hydrogen-like/helium- like 
line intensities for all ions measurable with the LETGS. 
As can be seen from Fig. the largest value for psi is 
obtained for Algol (psi = 1-54), while for the low-activity 
stars no lines from H-like or He-like silicon and magne- 
sium are observed. For Algol the p-values increase with 
decreasing atomic number reaching p^ — 8.37 for nitro- 
gen; no lines from carbon are observed. For Capella psi 
and pMg are below unity, then p increases to pn = 4.91 
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Fig. 1. Theoretical line fluxes for a cosmic abundance 
plasma for H-like and He-like resonance lines for C (solid 
lines), N (dotted lines), O (short-dashed lines), Ne (dash- 
dotted lines), Mg (dash-triple dotted lines), and Si (long- 
dashed lines) as calculated with the MEKAL code. Note, 
that for all elements the Ly a lines are stronger than the 
He-like resonance lines. 

for nitrogen while pc is lower. For the low-activity star 
Procyon all p- values except for carbon (pc = 1.78) are be- 
low unity; the neon lines are observed only quite weakly. 
So clearly the available data suggest large differences in 
the emission measure distribution for the sample stars. 

We will base our differential emission measure recon- 
struction solely on those lines. The disadvantage of us- 
ing those lines, i.e., their formation over relatively broad 
temperature range, is in our opinion more than compen- 
sated by a number of advantages: First of all, the atomic 
physics of hydrogen- and helium- like ions is much simpler 
than that of more complicated ions. Second, these lines are 
among the strongest lines; they can therefore be measured 
in a large sample of stars (cf. Ness et al. |2002b) and the 
DEM reconstructions of different stars can be compared 
with each other since they are computed in the same fash- 
ion. And third, as we will show below, these lines are very 
likely to already contain most of the temperature infor- 
mation contained in stellar coronae. It is important to re- 
alize that more line ratios do not necessarily provide more 
information on the temperature structure; additional line 
ratios may either contain no or no new information or may 
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Fig. 2. Theoretical line flux ratios of H-like by He-like 
resonance lines for C, N, O, Ne, Mg, and Si in comparison 
with the measured ratios for Algol, Capella, and Procyon. 

provide conflicting information. For example, Mcintosh 
( 200(J) gives the measured and calculated line ratios used 
in his differential emission reconstruction (cf . his Fig. 3b) , 
which deviate by almost an order of magnitude in the 
worst case. 

4. Which continuum to use ? 

In order to fix the overall normalization we use both the 
measured shape of the continuum and the absolute level 
of the observed continuum radiation. The first problem to 
be solved - a problem very familiar to optical astronomers 

- is the correct placement of the continuum. While strong 
lines can be recognized easily, the sum of weak lines, each 
of which remains undetected, can in principle produce a 
"pseudo-continuum". Since specifically the LETGS cov- 
ers such a large band pass it appears unlikely that over 
the whole instrument band pass from 5 A - 170 A such a 
"pseudo-continuum" is produced, while in narrower spec- 
tral bands this may well be the case. In order to isolate 
the continuum we use a median filter in the following way: 
In a predefined wavelength region - typically we use 0.5A 

- we calculate the median and use this value as the char- 
acteristic continuum level at that particular wavelength. 
Clearly, if too many lines are located in the wavelength bin 
considered, the thus derived continuum level is too high. 
This is specifically the case in the rather crowded region 
between 9 A and 18 A, where it is next to impossible to 
reliably place any continuum. Fortunately, other spectral 
regions are far less crowded and do allow a rather reliable 
continuum placement. 

5. A worked out example: Algol 

An 80 ksec observation of the eclipsing binary Algol has 
been carried out with the LETGS on board Chandra; the 
recorded data set and an analysis of the He-like and H- 
like lines has been presented by Ness et al. (|2()()2a|l . while 
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Schmitt and Ness I20Q2|) discuss the carbon and nitrogen 
abundances of Algol and other giants. Here we focus on 
the determination of the differential emission measure dis- 
tribution and elemental abundances; for a comparison of 
the coronal spectra of HR1099 and Algol B we refer to 
Drake <|2UtI3|l . 

5.1. DEM modeling with Chebyshev Polynomials 
5.1.1. Temperature structure 

Algol's X-ray emission is very strong and except for car- 
bon all H-like and He-like lines from N, O, Ne, Mg, and 
Si were detected; the failure to detect carbon lines in the 
X-ray spectrum of Algol is model-independent and due to 
a nitrogen enrichment of CNO-cycle processed material 
(Schmitt & Ness 120(13) . Thus unfortunately no informa- 
tion is available in the lower temperature range of the 
emission measure distribution from hydrogen- or helium- 
like lines. The flux ratio between the Si Ly a and He-like 
resonance line exceeds unity indicating that the peak of 
the emission measure distribution is beyond 10 MK. We 
therefore consider four values of r T maxl i.e., 20 MK, 30 MK, 
40 MK, and 50 MK. We first considered the simplest pos- 
sible combination of Chebyshev polynomials with M=4, 
constrained to yielding a positive emission measure dis- 
tribution. In order to prevent negative emission measures 
we introduced a penalty function resulting in large val- 
ues whenever the reconstructed emission measure distri- 
bution has negative values in the interval between T.,„j n 
and T m ax', we chose T m i n = 4 x 10 5 K, and note that 
our results are not very sensitive to the specific choice 
of T m i n . We specifically point out that the line with the 
coolest peak formation temperature is the resonance line 
of Nvi with a peak formation temperature of ps 10 61 K 
(cf. Fig. Q and with our Chandra LETGS data we have 
little information on the emission measure distribution be- 
low 10 6 . For each permitted emission measure distribution 
the resulting line ratios of N, O, Ne, Mg, and Si were com- 
puted and compared to the observed line ratios via the % 2 
test statistics. 

In this paper we use the Chianti software package (cf. 
Young et al. I2003[) to compute both line and continuum 
intensities for all modelling of Chandra LETGS data. The 
line intensities were computed in photon mode, for the 
continuum the contributions from bremsstrahlung, recom- 
bination continuum, and two photon continuum were sep- 
arately computed and added. The ionization equilibrium 
by Mazzotta et al. I|1998fl was used. All calculations were 
carried out using cosmic abundance as quoted by Allen 
l|1973J) . For the relevant elements we specifically used the 
values [C/H]= 8.52, [N/H]= 7.96, [0/H]= 8.82, [Ne/H] = 
7.92, [Mg/H]= 7.42, [Si/H]= 7.52, and [Fe/H] = 7.60. 
Comparing these abundance values to the most recent val- 
ues given by Grevesse and collaborators (1998), we find 
essentially identical values for C, N, O and Si, while for 
the elements Ne, Mg and Fe Grevesse & Sauval Ijl998[l 
give [Ne/H]= 8.08, [Mg/H]= 7.58 and [Fe/H] = 7.50, i.e., 
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Fig. 3. Best fit differential emission measure distributions 
with fourth-order Chebyshev polynomials and tempera- 
tures T max = 20 MK, 30 MK, 40 MK, and 50 MK; only 
the choices T max = 30 MK and T rnax = 40 MK yield ac- 
ceptable fits to the Chandra LETGS data. 

values differing by factors of 1.44 (Ne and Mg) and 0.79 
(Fe) respectively. We only give the multiplicative factors, 
by which our derived abundance values differ from the 
ones quoted by Allen 1|1973[) : in order to convert to the 
ones quoted by Grevesse & Sauval ( 1998), the He and Mg 
abundances must by multiplied by 0.69, the Fe abundance 
by 1.26. 

Our modeling attempts showed that already with the 
choice of M = 4 good fits to the line ratios of the Ly Q and 
He- like resonance lines can be obtained. In Fig. [3] we plot 
the best fit reconstructed emission measures (for the case 
M = 4) for the peak temperatures T max = 20 MK, 30 MK, 
40 MK, and 50 MK; both the curves for T max = 30 MK 
and 40 MK yield acceptable fits (solid lines), while the 
choices of T max = 20 MK and T max = 50 MK lead to 
unacceptable fits (cf. Tab. [I}. Our formal fit results are 
presented in Tab. ^ where we give for those four best 
fits the resulting values of x 2 f° r the line ratio fit as well 
as the modelled line ratios. A maximum temperature of 
Tmax = 20 MK is simply too low to explain the observed 
ratio between Ly Q and He-like resonance line for silicon. 
On the other hand, for the model with T max — 50 MK 
and M=4 too much emission measure is located at high 
temperatures leading to an increase in \ 2 - I n au fits we 
also checked for the goodness of fit to the continuum; in 
general, fits with higher temperature result in better con- 
tinuum fits than lower temperatures. 

Clearly, the restriction to the first four Chebyshev 
polynomials results in rather simple emission measure dis- 
tributions, which, however, provide good line ratio fits 
for the correct choices of T max . How unique are the de- 
rived emission measure distributions ? In order to assess 
this issue we introduced higher order Chebyshev polyno- 
mials, which result in more complicated emission mea- 
sure distributions and in improved fits. In Fig. 0] we plot 
the best fit reconstructed emission measures for the cases 
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Table 1. Line ratio fit results for DEM models based on Chebyshev polynomials. The numbers given 
in parentheses are (T max , M). For all analysed atomic species we give the observed line ratios with their errors as well 
as the modelled line ratios for the various models considered. The last row gives the values of the x 2 test statistics for 
the respective models. 



Element 


Observed 


Modelled 


Modelled 


Modelled 


Modelled 


Modelled 




ratio/error 


line ratio 


line ratio 


line ratio 


line ratio 


line ratio 






(20,4) 


(30,4) 


(40,4) 


(50,7) 


(50,4) 


N 


8.38±1.28 


8.85 


9.44 


9.21 


9.00 


6.66 





8.00±0.71 


7.37 


7.02 


7.04 


6.84 


5.94 


Ne 


4.46±0.28 


4.80 


4.61 


4.76 


4.41 


4.77 


Mg 


2.41±0.31 


2.27 


2.52 


2.69 


2.48 


2.97 


Si 


1.55±0.14 


1.10 


1.49 


1.66 


1.63 


1.97 


x 1 




12.45 


3.16 


4.82 


3.36 


12.89 



/ : 



6 r 




Temperature (K) 



Fig. 4. Best fit differential emission measure distribu- 
tions with Chebyshev polynomials with (T max = 30 MK, 
M=4; dashed curve), {T max = 40 MK, M=4; solid curve), 
and (T max = 50 MK, M=7; dotted curve). All parameter 
choices yield acceptable fits to the Chandra LETGS line 
ratio data, the higher temperature yield better continuum 
values. 

{T max = 30MK,M = 4), (T max = 40MK,M = 4) and 
(T m ax — 50MK,M = 7) as a function of temperature; 
the resulting y 2 values are given in Tab. As can be 
seen from Tab. ^ au fits have similar goodness of fit pa- 
rameters, but the resulting DEM curves are quite dif- 
ferent. In particular, assuming T max — 50 MK, leads to 
a bimodal emission measure distribution with a second 
peak at 45 MK, corresponding to a cutoff energy of about 
3keV. The presence of such a peak in the emission mea- 
sure distribution can be readily recognized from the high 
energy continuum emitted by hot plasma. However, the 
LETGS is not particularly sensitive in that wavelength 
range. For XMM-Newton data, for example, simultane- 
ously taken CCD spectra at higher energies would yield 
important constraints at the high temperature end of the 
DEM distribution which are not provided by the LETGS. 
From a statistical point of view the improvements in fit 
quality are so small that they do not warrant the intro- 
duction of additional degrees of freedom. 
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Fig. 5. Chandra LETGS spectrum for Algol binned in 
wavelength bins with a width of TBD A(dots) and best fit 
differential emission measure model continuum for M = 4 
and Tmax = 40 MK (solid line); dashed curves indicate a 
15% systematic uncertainty in effective instrument areas. 
Note that the continuum has been calculated from the 
abundances actually derived for Algol. 

In order to demonstrate the overall goodness of fit 
of our continuum models we plot in Fig. [S] the observed 
Chandra LETGS for Algol (with a resolution of 0.03 A; 
small dots) and the best fit continuum model (solid line) 
for the case T ma x = 40 MK and M = 4. We empha- 
size that in fixing the normalization we attempted only 



10 



J.H.M.M. Schmitt and J.-U. Ness: Coronal abundances from high-resolution X-ray data: The case of Algol 



modeling the continuum, but not the individual lines; the 
continuum modeling includes higher orders up to order 
ten. The continuum was computed with the specific set 
of abundances calculated for Algol, but no iterative (re- 
determination of abundances was performed, since the 
changes in the continuum level are of the order of a few 
percent at most. The continuum was fitted in the range 
19 A - 40 A and 70 A - 105 A. Fig.0shows that the contin- 
uum is well fitted at short wavelengths and in the wave- 
length range 80 - 100 A. The fit is extremely poor in the 
wavelength range between 8 - 18 A; this is hardly surpris- 
ing since numerous emission lines are located (cf. Fig. 1 
in Ness et al. I2002afl in that wavelength region. Our fit 
describes quite well the carbon edge at 44 A. Our fits be- 
have badly near « 53 A and 63 A, where the LETGS spec- 
tra show two dips, which are instrumental and caused by 
gaps in the HRC-S microchannel plates. At long wave- 
lengths our continuum fit describes the observed data, but 
is somewhat high. We checked our analysis procedures on 
the public Chandra LETGS data (500 ksec) of the isolated 
neutron star RXJ1856. 5-3754 and found very good agree- 
ment between our results and those published in the liter- 
ature. We have at present no satisfactory explanation for 
the discrepancy at longer wavelengths. On the one hand 
our choice of median filtering has some bias towards higher 
values because the filtering procedure essentially regards 
spectral lines as background fluctuations. On the other 
hand errors in the instrument calibration can also not be 
ruled out; calibrations in the EUV are notoriously difficult 
and an error in the relative calibration between the long- 
and short wavelength region will also help to reduce the 
observed discrepancy. An even larger absorption column 
towards Algol would also improve the fit. 

5.1.2. Abundances 

With the reconstructed abundance independent emis- 
sion measure distributions (cf. Fig. |3J) and the measured 
line fluxes we can now calculate elemental abundances. 
Specifically, we can calculate for each observed spectral 
line that value of abundance which forces agreement be- 
tween predicted and observed line fluxes. Obviously, the 
abundances for a given chemical element will depend on 
the specific lines used and any errors in the atomic param- 
eters and temperature structure will propagate into the 
derived abundances; our results for various lines are sum- 
marized in Tab. [21 All abundance values quoted are rela- 
tive to the values [C/H]= 8.52, [N/H]= 7.96, [0/H]= 8.82, 
[Ne/H]= 7.92, [Mg/H]= 7.42, [Si/H]= 7.52, and [Fe/H] = 
7.60. An inspection of Tab. [3 shows that H- and He-like 
lines give consistent abundance values; this is of course 
no surprise, since the adopted differential emission mea- 
sure distribution fits the line ratios very well. The upper 
limit derived from the absence of the C vi Ly a line im- 
plies that the carbon abundance (relative to solar) is at 
the level of less than 10 percent, i.e., most of the car- 
bon has been depleted. No iron lines have been used for 



the DEM reconstruction. Inspection of Tab. shows that 
the nitrogen abundances for accpetable models lie be- 
tween 1.9-2.9, the oxygen abundances between 0.25-0.36, 
the neon abundances between 0.95-1.34, the magnesium 
abundances between 0.48-0.68, and those of silicon be- 
tween 0.42-0.61. The Fe abundances (for the best fit DEM 
models) vary from 0.10 - 0.49 depending on the lines cho- 
sen. These discrepancies are of course worrisome, but they 
are small compared to the discrepancies found in solar 
work (cf. Mcintosh, 2000). Also, some Fe lines consistently 
yield abundance values inconsistent with those derived 
from other lines. For example, the FeXVll 17.07 A line 
yields larger values than the neighboring 15 A line, and 
the Fexix 101.55 A line gives much higher values than 
the other Fe XUV lines, indicating systematic problems 
with those lines. Tab. also reveals that the Fe abun- 
dances derived from XUV lines (with the exception of 
Fexix 101.55 A)are in general lower than those of X-ray 
lines. Whether this is due to atomic physics or instrumen- 
tal problems is difficult to say at present. The errors in 
the effective area calibration of the LETGS are thought 
to be of the order ± 15 % at most; the magnitude of this 
effect appears to be too small to account for the observed 
discrepancies. Another issue is the value of the interstellar 
absorption column density; for our analysis we used the 
value of 5xl0~ 18 cm~ 3 , and higher values would reduce, 
smaller values increase the discrepancy between model and 
observations. 

5.2. Temperature component analysis 
5.2.1. Temperature structure 

It is of course also possible to carry out the analysis of the 
line ratio data in a more traditional way using individ- 
ual, discrete temperature components. In this approach 
one uses L independent spectral components with tem- 
peratures 7], I = 1 . . .L and emission measures A\. One 
clearly ought to demand 

Ai>0, (32) 

and, since we are using line ratios, any solution is only 
determined up to a multiplicative factor which needs to 
be fixed by the continuum. We use the constraint 

Ef =1 ^ = l (33) 
and therefore a solution with L spectral components has 
2L — 1 fit parameters. In Tab.^lwe present the results of 
our analysis of the available line-ratios for Algol using in- 
dividual temperature components; we list the best fit tem- 
peratures and amplitudes as well as the goodness of fit pa- 
rameter x 2 ■ An inspection of Tab. shows that a good fit 
to the observed line ratios measured for the Ly a - and He- 
like resonance lines for N, O, Ne, Mg, and Si requires four 
spectral components with a total of seven adjustable pa- 
rameters. Descriptions with fewer components and fewer 
adjustable parameters cannot explain the observed line ra- 
tios; the test statistic \ 2 depends - not surprisingly - sen- 
sitively on the number of spectral components, but three 
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Table 2. Abundance determinations for DEM models based on Chebyshev polynomials; listed are the elements and 
lines used (first two columns), the measured number of counts, the modelled number of counts as well as the abundance 
relative to the Allen l|1973|l abundance values for the various models characterized by number of polynomials and peak 
temperature. Last row indicates whether fits are acceptable "acc" or unacceptable "not acc" . 



Ion 


Line 

(A) 


Obs. 


Mod. 


(7,50) 


Mod. 


(4,30) 


Mod. 


(4,40) 


Mod. 


(4,50) 


Mod. 


(4,20) 


Cvi 


33.74 


< 60 


588.2 


< 0.10 


763.1 


< 0.08 


660.5 


< 0.09 


659.8 


< 0.09 


908.1 


< 0.05 


Nvi 


24.785 


1119.1 


416.9 


2.68 


567.4 


1.97 


501.1 


2.23 


467.7 


2.39 


673.0 


1.71 


Nvn 


28.79 


141.3 


48.6 


2.91 


63.3 


2.23 


57.0 


2.48 


54.4 


2.59 


76.9 


1.84 


vin 


18.97 


2883.0 


7925.2 


0.36 


10289.0 


0.28 


9057.4 


0.32 


8890.1 


0.32 


11749.0 


0.25 


Ovn 


21.6 


262.5 


836.7 


0.31 


1059.7 


0.25 


930.1 


0.28 


938.5 


0.28 


1251.5 


0.21 


Nex 


12.138 


2481.5 


1956.4 


1.27 


2602.2 


0.95 


2300.1 


1.09 


2194.7 


1.13 


2996.8 


0.83 


Neix 


13.45 


665.0 


497.8 


1.34 


633.5 


1.05 


542.4 


1.23 


558.4 


1.19 


710.7 


0.94 


Mg xii 


8.43 


578.0 


878.0 


0.66 


1196.3 


0.48 


1086.1 


0.53 


985.0 


0.59 


1406.1 


0.41 


Mgxi 


9.17 


224.1 


329.3 


0.68 


439.7 


0.51 


374.9 


0.60 


369.4 


0.61 


514.6 


0.44 


Si xiv 


6.15 


658.3 


1135.8 


0.58 


1504.8 


0.44 


1452.2 


0.45 


1274.1 


0.52 


1784.5 


0.37 


Si xin 


6.65 


480.7 


784.7 


0.61 


1132.1 


0.42 


983.0 


0.49 


880.3 


0.55 


1618.6 


0.41 


Fe xvn 


15.01 


1018.4 


5322.6 


0.19 


6740.0 


0.15 


5714.8 


0.18 


5970.6 


0.17 


7286.7 


0.14 


Fexvn 


15.27 


364.7 


1286.1 


0.28 


1628.1 


0.22 


1380.6 


0.26 


1442.6 


0.25 


1760.1 


0.21 


Fexvn 


17.07 


1124.5 


2278.8 


0.49 


2879.7 


0.39 


2443.5 


0.46 


2556.2 


0.44 


3109.9 


0.36 


Fexvn 


93.92 


258.0 


2023.9 


0.13 


2599.7 


0.10 


2196.2 


0.12 


2270.3 


0.11 


2857.9 


0.09 


Fexix 


108.37 


203.5 


1193.6 


0.17 


1571.0 


0.13 


1323.1 


0.15 


1339.0 


0.15 


1817.1 


0.11 


Fexix 


101.55 


180.9 


444.9 


0.41 


585.5 


0.31 


493.2 


0.37 


499.0 


0.36 


677.3 


0.25 


Fexxi 


128.73 


266.32 


1914.5 


0.14 


2725.2 


0.10 


2300.5 


0.12 


2147.8 


0.12 


3730.5 


0.07 


Fe xxiii 


132.85 


1178.6 


5561.9 


0.21 


8505.0 


0.14 


7508.5 


0.16 


6250.8 


0.19 


11021.3 


0.11 






acc 


acc 


acc 


not acc 


not acc 



components, yielding a test statistic Xmin = 17-15, are 
just not sufficient, while the introduction of a fourth and 
fifth component reduces the level of x 2 1° below unity. 
The necessity to introduce 7 parameters in order to fit 5 
data points is a nightmare for every statistician and one 
wonders about the effective number of degrees of freedom. 
Presumably these problems have to do with the positiv- 
ity constraint Eq. |23 which severely limits the available 
solution space. Unfortunately, the resulting best fit tem- 
perature quadruple is not unique. If one accepts a limit 
of xtim < 6 as a threshold below which solutions are ac- 
cepted uncertainties of up to 0.2 dex in Log T result. 

5.2.2. Abundances 

One expects that the derived abundances depend sensi- 
tively on the choice of the adopted temperature compo- 
nents if one is working with a small number of temperature 
components. These expectations are verified by an inspec- 
tion of Tab. |31 where we list the abundances derived from 
our temperature fits. In Tab. [3] we list the derived elemen- 
tal abundances for each discrete temperature solution for 
three, four, and five individual components (list in the first 
row). The depletion of carbon is recognized independent 
of the chosen number of temperature components, while 



abundances for individual elements and lines can but must 
not vary considerably. For example, the O abundance de- 
rived from Ovni 18.97 A changes from 0.19 - 0.20, i.e., is 
essentially model independent, while the iron abundance 
derived from Fe xx 101.55A varies from 0.06 to 0.24 in an 
extremely model dependent way. 

If we now compare the abundances from the simplest 
acceptable discrete temperature component model (i.e., 
4-T) to those derived from the simplest continuous emis- 
sion measure distribution model (4, 40) we find that the 
nitrogen and oxygen abundances differ by a factor of ss 
2, the discrete temperature component abundances being 
lower. For the neon, magnesium, and silicon abundances 
the situation appears similar with the continuous temper- 
ature distribution abundances being higher, while for iron 
the two sets of abundances agree reasonably well the main 
difference coming from the lines used for abundance anal- 
ysis. 

5.3. Towards a physical model 

It should have become clear that neither the modeling ap- 
proach based on Chebyshev polynomials nor that based 
on individual temperature components includes a great 
deal of physics other than the theory of hot thermal 
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Table 3. Abundance determinations for discrete temper- 
ature component models; listed are the elements and lines 
used (first two columns), the model parameters in paren- 
theses indicating number of temperature components and 
value of x 2 as weu as the individually derived abundances 
relative to cosmic abundances. Last row indicates whether 
model gives acceptable "acc" or unacceptable fits "not 
acc" . 



Ion 


Line (A) 


(3,9.1) 


(4,2.26) 


(5,0.83) 


C vi 


33.74 


< 0.06 


< 0.06 


< 0.06 


N vi 


28.79 


1.54 


1.41 


1.45 


N vii 


24.785 


1.35 


1.41 


1.45 


O vii 


21.6 


0.19 


0.17 


0.19 


O vm 


18.97 


0.20 


0.19 


0.19 


Neix 


13.45 


0.68 


0.77 


0.74 


Ncx 


12.138 


0.58 


0.69 


0.69 


Mgxi 


9.17 


0.28 


0.33 


0.33 


Mgxn 


8.43 


0.39 


0.36 


0.37 


Si xm 


6.65 


0.33 


0.35 


0.36 


Si xiv 


6.15 


0.34 


0.34 


0.34 


Fe xvn 


15.01 


0.06 


0.09 


0.09 


Fe xvn 


15.27 


0.09 


0.13 


0.13 


Fe xvn 


17.07 


0.16 


0.24 


0.23 


Fe xvn 


93.92 


0.05 


0.05 


0.05 


Fexx 


101.55 


0.24 


0.15 


0.06 


Fexix 


108.37 


0.10 


0.06 


0.06 


Fexxi 


128.73 


0.16 


0.14 


0.13 


Fe xxiii 


132.85 


0.11 


0.15 


0.15 






not acc 


acc 


acc 



The parameters a, (3, and 7 in Eg. 1341 are the power law 
coefficients in the laws for radiative cooling, loop heating, 
and loop cross section expansion which are assumed to 
have the form 



A(T) = A T 7 Q(T) = Q T> 3 A(T) = A T a 



(35) 



respectively. Clearly, these assumptions are somewhat ide- 
alistic, and in particular the assumption about the heating 
law is totally arbitrary, but the resulting structure is found 
to be remarkably insensitive to the detailed form of these 
coefficients. The constant C norm in Ea . 1341 is given by the 
expression 



'kq (2a + 7 + 3/2) 



8A fc 2 



(36) 



where kq denotes the constant in the heat conduction law. 
The singularity at the loop top, where T = T max , arises 
from the boundary condition = at the same place. 
In such structures the fundamental physical parameters 
pressure, loop top temperature, and the loop length L are 
not independent but are related through so-called scaling 
laws viz. 



p L ~ T„ 



(37) 



In order to make progress we assume ad hoc that the dis- 
tribution of the actual loop top temperatures $(T m(II ) of 
a collection of loops can be described through a Gaussian 
viz. 



plasma emission. Let us consider a "standard" solar loop 
with apex temperature T max = 2 x 10 6 K, pressure p = 
ldyncm -2 , and loop half length L = 2 x 10 9 cm. Such 
a structure contains a total emission measure of a few 
times 10 46 cm -3 . Comparing this estimate to the total so- 
lar coronal emission measure or to the emission measure 
observed from stars it is clear that several hundreds and 
possibly more than thousand such " standard" loops must 
be contributing to the emission observed at any given 
time. This is amply demonstrated by the thousands of 
YOHKOH and SO HO images, which show a vast variety 
of different emission structures, and only during a stronger 
flare an individual structure can dominate the total X- 
ray output. In a stellar context the situation is less clear. 
Again, during a large flare the overall emission is certainly 
dominated by the emission from the flare region alone, for 
the quiescent emission we assume that a larger number of 
individual features is responsible for the observed emis- 
sion. 

If one assumes that the X-ray emission comes from 
a sample of individual constant pressure "atmospheres" 
extending from some lower temperature Tq to some maxi- 
mum temperature T max , the differential emission measure 
distribution £(T) of such a loop can be calculated as (cf. 
Bray QHD) 



${T max ) 



(38) 



with a mean loop top temperature Tq and some dispersion 
ot ■ We further assume - for the time being - that all loops 
have the same length. It is clear that only loops with loop 
top temperature T max > T can contribute to the overall 
differential emission measure at some given temperature 
T. The total differential emission measure from such a 
collection of loops is therefore given by 



max- 1 - max 



2- 7 +/3 



(39) 



where we have used the scaling law Eq. 1371 to replace the 
(constant) loop pressure p with T max and L. Because of 
the scaling law Ea. 1371 v and T max are not independent, 
and in fact for given L the pressure and thus the total 
emission measure depend quite sensitively on T max . If we 
then - finally - assume that the loop length L is related 
to the loop top temperature with some power law coeffi- 
cient (5, we are led to the following analytical form of the 
differential emission measure distribution: 

Ctot\-F) — A norrn T 
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with the three power law coefficients a\, a%, and 023, the 
mean loop top temperature T and its variance ot as pa- 
rameters describing the shape of the distribution function. 
With the Chandra LETGS data we can test whether the 
shape specified by Eq. 0U] is consistent with our observa- 
tions. Calculations showed that the hydrogen- to helium- 
like line ratios depend sensitively on a±, Tq, and <jt, but 
less sensitively on the parameters a 2 and 013. The best 
fit model was found to be given by the set ai = 0.48, 
a 2 = 0.50, a 3 = 0.50, T = 7MK, and a T = 11 MK. In 
Fig. we plot the DEM distribution function in compari- 
son to the best-fit DEM distributions reconstructed from 
fourth-order Chebyshev polynomials with T max — 30 MK 
and T max = 40 MK; as is clear from Fig. El the DEM 
reconstructions agree extremely well with each other for 
temperatures above 1 MK, while the slope of the emission 
reconstructions disagrees for lower temperatures. Since 
there are almost no constraints on the DEM in this tem- 
perature range (note that the peak formation temperature 
for N VII is at 1 MK), such a discrepancy is hardly surpris- 
ing. As a first conclusion we therefore note, that a function 
form as given in Eq. 1401 provides a reasonable description 
of the observed line ratios with suitably chosen parameters 
To, <tt, cti, 0:2, and 03. In the framework of the chosen 
model it is a little surprising that the mean (Gaussian) 
loop top temperature is rather low (T = 7MK). However, 
because of the large dispersion of 11 MK, significantly 
hotter loops occur with high probability. Since accord- 
ing to the scaling laws higher loop top temperatures im- 
ply higher pressure, higher pressure implies higher density, 
and higher density implies larger emission measure, so the 
hotter loop will contribute most. For constant cross sec- 
tion loops with the canonical radiative loss function one 
expects ai ~ 0; the observed value a\ ~ 0.5 can only be 
explained by expanding loop geometries. 

6. Discussion and Conclusions 

The new generation of X-ray spectrometers on board the 
Chandra and XMM-Newton satellites allows the deter- 
mination of elemental abundances in hot X-ray emitting 
plasmas. The Chandra LETGS has the specific advan- 
tage of a very large band pass with an ensuing sensi- 
tivity to lines from rather different temperature regimes. 
Our analysis of the Algol LETGS spectrum shows that 
the abundances for the elements neon, magnesium, sili- 
con, and iron are all sub-solar. This is in line with previ- 
ously published high-resolution abundance determinations 
of HR1099 with the Chandra HETGS (Drake et al- HSD 
and XMM-Newton (Brinkmann et al. l2001fl . Among those 
elements neon has the relatively highest abundance, i.e., it 
is least sub-cosmic. It appears that these conclusions are 
rather robust, and specifically do not sensitively depend on 
the methodology used ("global fit" vs. "emission line anal- 
ysis"). However, an inspection of Tab. [2] shows that, for 
example, the iron abundance determinations considerably 
depend on the lines used for the analysis. The short wave- 
length lines of Fexvn at 15.01 A, 15.27 A, and 17.07A 




10 5 10 6 10 7 



Temperature (K) 

Fig. 6. Best fit differential emission measure distribution 
derived from Ea.l40lfsolid curve) compared to DEM distri- 
butions derived from fourth order Chebyshev polynomials 
with T max =30 MK and T max =40 MK. 



yield higher abundances than the Fe lines at 93.92 A, 
108.37 A, 128. 37 A, and 132.82 A; the Fe line at 101.55 A 
yields larger abundances than the rest of the EUV Fe 
lines. The reason(s) for this discrepancy are not quite 
clear. Optical depth effects in the Fexvn are very likely 
not the cause (Ness et al. 2003J) and would in fact even 
worsen the discrepancy. Since the short wavelength lines 
are all from Fexvn, the calculated emission measure dis- 
tribution for this ion might be incorrect, i.e., too large. 
Since Fe xvn is produced over a rather large temperature 
range this appears somewhat unlikely since then also the 
continuum emission would have to be incorrectly placed. 
The Fe xvn long wavelength lines are affected by absorp- 
tion, but we used already a rather large value of Njj', low- 
ering the absorption column density would again worsen 
the discrepancy. Systematic errors in the instrument cal- 
ibration might affect the long wavelength portion of the 
spectrum in a different way than the short wavelength 
portion, but the magnitude of the effect is much larger 
than the systematic calibration uncertainties (~ 15 %). 
Finally, atomic physics uncertainties might affect the long 
wavelength lines different compared to the short wave- 
length lines. At any rate, we have to conclude and state 
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that we presently have no satisfactory explanation for the 
discrepant abundance determinations for iron. 

As to the abundance discrepancies derived from dif- 
ferent analysis methods, our studies have clearly demon- 
strated the importance of the correct determination of the 
underlying temperature structure for a correct determina- 
tion of elemental abundances. The cooling functions of in- 
dividual lines contribute significantly over a temperature 
range of 0.3 dex and the shape of the emission measure dis- 
tribution also implies considerable contributions far away 
from the peak formation temperatures of individual lines. 
The lines used in our study, i.e., Ly Q and He-like r lines, 
are formed over a rather wide temperature range, other 
lines, in particular lines from ions with incompletely filled 
shells, are formed over somewhat narrower temperature 
ranges necessitating an even better knowledge of the tem- 
perature structure. We purposely used only those lines for 
our differential emission measure reconstruction, because, 
first, these lines are among the strongest for each element 
and therefore the most likely lines to be detected in a 
recorded X-ray spectrum, and second, the atomic physics 
of those lines ought to be known best. A reliable method 
for abundance determination must prevent any cross talk 
between the temperature and abundance structure of a 
plasma; therefore, the temperature structure should be 
determined independent from the elemental abundances 
either from line ratios of lines of the same chemical ele- 
ment (as done in this paper) or by using lines only from 
the same element (as done by, e.g., Drake et al. I2001|l . 
We next emphasize the need of physical considerations in 
the determination of the temperature structure. This is 
in particular required if one ever desires to determine el- 
emental abundances in the X-ray range with an accuracy 
achieved by optical abundance determinations. We have 
at the moment only few clues as to the differential emis- 
sion measure distributions realized by stellar coronae and 
the uncertainties in our knowledge of the correct temper- 
ature structure prevents from reaching precisely this goal. 
A modeling of the coronal emission in terms of individual 
temperature component is unsatisfactory from a physical 
point of view, from a procedural point of view and from 
a mathematical point of view. Abundances determined in 
this way may have small statistical errors (of a few per- 
cent depending on the SNR of the modeled data), but 
rather large systematic errors of 100 % or more; neverthe- 
less they are adequate to reveal general trends in abun- 
dance patterns. This is exemplified in Tab. 01 which shows 
that for example the oxygen abundance changes by a fac- 
tor two for models with discrete temperature components. 
A comparison of the abundances derived for Algol in this 
paper with those derived by Antunes et al. I| 1994(1 from 
ASCA using a two-temperature variable abundance mod- 
eling approach also shows that the general trend in the run 
of elemental abundances is captured and the "low" iron 
abundance and "high" neon abundance are recognized, 
while the abundances of individual elements can vary by 
at least a factor of two. Also, the real clue of the Chandra 
LETGS Algol observation, i.e., the overabundance of ni- 



trogen with its profound physical implications (cf. Schmitt 
and Ness 12002|l went unnoticed in the modeling with the 
lower resolution ASCA data. 

As to Algol specifically, our detailed temperature and 
abundance modeling confirms the results previously de- 
rived by Schmitt and Ness (|2~0~02"|1 . Because of temperature 
dependence of the emissivity functions of the Ly a -lines for 
C and N (cf. Fig. QJ, the line ratio between these lines 
must stay below 0.57 (for cosmic abundances) regardless 
of the underlying temperature structure in contrast to the 
observed ratio of > 23. Our modeling now shows that car- 
bon is depleted down to at least 8%, while nitrogen is 
enhanced by about 70% or more (all relative to cosmic 
abundances). This effect is dramatic. Assuming the cos- 
mic abundance pattern recommended by Holweger (2001 ) 
there are 4.58 carbon atoms for every nitrogen atom, while 
in Algol's corona we have (at least) 8.2 nitrogen atoms for 
every carbon atom ! This reversal of carbon and nitrogen 
abundance can be readily explained by assuming that one 
is studying CNO-cycle processed material in the corona of 
Algol B, since the equilibrium abundance of CNO nuclei 
participating in the cycle is such that most nuclei occur 
as N 1A nuclei (Caughlan 1965). In no other spectral range 
than the X-ray band can the chemical abundance of the 
B component of the Algol system be studied. 

Our Chandra LETGS spectrum of Algol thus demon- 
strates the wealth of physical information contained in an 
X-ray spectrum with high spectral resolution and - at the 
same time - good signal-to-noise ratio. The latter is as im- 
portant as the former, since data with poor signal-to-noise 
will not allow the derivation of meaningful and significant 
results. The exposure of such spectra requires substantial 
satellite resources, yet it represents the only way to extract 
information on the physics of stellar coronae. 
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High State 
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